


Institutional Archive of the Naval Postgraduate School 





Calhoun: The NPS Institutional Archive 
DSpace Repository 


Theses and Dissertations 1. Thesis and Dissertation Collection, all items 


1969-12 


An experimental single Fourier transform 
optical processor applied to film grain analysis 
and a ship wave detection method 


Ferguson, Robert Dale 


Monterey, California. U.S. Naval Postgraduate School 
http://ndl.handle.net/10945/12177 


This publication is a work of the U.S. Government as defined in Title 17, United 
States Code, Section 101. Copyright protection is not available for this work in the 
United States. 


Downloaded from NPS Archive: Calhoun 


Calhoun is the Naval Postgraduate School's public access digital repository for 
(8 DUDLEY research materials and institutional publications created by the NPS community. 
«ist sae Calhoun is named for Professor of Mathematics Guy K. Calhoun, NPS'‘s first 


INN KNOX appointed — and published -- scholarly author. 

| LIBRARY Dudley Knox Library / Naval Postgraduate School 

411 Dyer Road / 1 University Circle 
Monterey, California USA 93943 





http://www.nps.edu/library 





NPS ARCHIVE 
1969 
FERGUSON, R. oo : 


AN EXPERIMENTAL SINGLE FOURIER TRANS- 

FORM OPTICAL PROCESSOR APPLIED TO FILM 

GRAIN ANALYSIS AND A SHIP WAVE DETEC- 
TION METHOD 


by 


Robert Dale Ferguson 





DUDLEY KNOX LIBRARY 
NAVAL POSTGRADUATE SCHOOL 
MONTEREY, CA 93943-5101 


United States 
Naval Postgraduate School 


AN EXPERIMENTAL SINGLE FOURIER TRANSFORM 
OPTICAL PROCESSOR APPLIED TO FILM GRAIN 
ANALYSIS AND A SHIP WAVE DETECTION METHOD 


Robert Dale Ferguson 





December 1969 


This document has been approved for public re- 
Lease and sake; its distribution 46 unloncted. 


a £ ¢*. ak ’ Ps 
; J , é 2&8, 


LIBRARY ” 
wiVAL POSTGRADUATE SCHOOL 
pwNTLNIY, CALIF. 98940 








An Experimental Single Fourier Transform Optical Processor 
Applied to Film Grain Analysis and A Ship Wave Detection Method 


by 


Robert Dale Ferguson 
Lieutenant, United States Navy 
B. S., University of Washington, 1963 


Submitted in partial fulfillment of the 
requirements for the degrees of 


ELECTRICAL ENGINEER 
and 
MAS TE RoC) eS Gul C bel Eb ip amuriGrAie: NGL NER RING 


From the 


NAVAL POSTGRADUATE SCHOOL 
December 1969 


ABSTRACT 


A single Fourier transform optical system was constructed and tested 
on two projects. First, film grain in a single Polar transform operation 
was studied and a mathematical model of the film grain and a semi-quan- 
titative prediction of the grain's transform were developed and verified 
using the experimental system. The film grain was found to have a noise- 
like spatial frequency spectrum which — a decreasing function of fre- 
quency, end that this spectrum could be considered Gaussian only in 
the case of fine grain films containing binary signal images. Second, 
the system was uSed to investigate.a proposed method of detecting a 
ship wave. A deterministic relation between a ship wave and the Fourier 
transform of its photograph was shown to exist. A Sample photograph of 
a ship wave was inserted into the system and operated on by a matched 
filter rotated in a plane normal to the optical axis, The output light 
intensity was determined as a function of the filter angular position. 

The output intensity when considered to be periodic was found to have a 


discernible spectral content due to the ship wave. 
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I, INTRODUCTION 


A, BACKGROUND 

Optical systems have received increased attention from physicists 
and engineers. The reason lies in the Fourier transform properties of a 
coherent optical system, and the laser has recently provided a practical 
source of coherent light. The optical image possesses two degrees of free- 
dom (i.e.; the coordinates of a point on the surface) in contrast to the 
Single degree of freedom, time, as found in electrical systems. With 
two degrees of freedom, an optical system appears to offer distinct advan- 
tages over an equivalent electrical system in some applications. 

Since the electrical engineer is generally not familar with the 
Fourier properties of an optical lens, a brief description is presented 
here. More detailed analyses are given in Preston [1], Francon [2], 
and Goodman [3]. Consider the system diagrammed in Figure 1. The 
input plane, Pl, is illuminated by a planar, unpolarized, unmodulated, 
monochromatic light wave described by the equation: 
f,(t) = f.e7%" (1-1) 
Where w is the radian frequency of the light wave. The plane, Pl, 
contains a light modulator t(x,,y,) such as a photographic transparency, 
The transmitted light is: 

i, yy a fae (x, ,y,)e* (1-2) 

where t is complex, Equation (1-2) may be written as: 


ims, Yuu = f t(x,,y,J)expljat+ é ev. )8 (1-3) 
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Assume the lens to be ideal and the aperture, lavge so that vignetting 
does not occur. Geometric optics shows that the energy diffracted at Pl 
in the direction, v, will be summed at a point in the back focal plane, P2. 
It may be shown that the phase at the summing point, (Xz,yz), is: 

¢ (Xo ,Ya) = d(x,y,)+ + rr (x. Xgty, Yo) — (1-4) 
where Ais the wavelength of light, Ris shown in Figure 1, and dis the 
optical path length from the origin of Pl to the point (xz,y,). The light 


incident at the point, (Xz,ys), becomes: 


j ast . 
iy bey Yo: t) =f’ ff t (x, t y,)exp J L g(x, ‘ i), Ue 
Pi 


a eso) 7 Cj} dx,dy, (1-5) 


The phase modulation term, ( (x,,¥,), is usually made constant in prac- 
tice and most often detection is by energy so the term outside the integral 
may be replaced by its rms value. For small diffraction angles where R 

is approximately equal to the focal length, F, equation (1-5) reduces to: 


Re 4 
fo(x5, yz) =f, ii tGowaypexp j \F (X,Xot yi V2)dx,,dy, (1-6) 
ea 


where f =Feexp 1/4 ( X1,¥,)+ — 
Equation (1-6) is the two-dimensional, bounded, Fourier transform of the 
function, t(x,y,), multiplied by a system constant. This may be more 
readily apparent by rewriting the kernel function exp j(27/A F) (x, Xsty,ye) 


as exp -j (WX +WV ) and defining: 


~aae 27 
ee oT, “ 
ast aa SUE Yo (1 7) 
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Thus Qi. and uw). are spacial frequencies defined in the back focal plane, 


Y 
P2. Two successive optical transform operations will map the input image 


into the output with only a coordinate reversal and a constant amplitude 


change. 


B. APPLICATIONS 

One of the more obvious applications follows directly from the above. 
Since the plane of the first transform, or frequency, plane displays the 
spatial frequency content of the input function, t(x,,y,), frequency 
modifying operations such as band limiting and matched filtering may be 
performed in this plane. The theory is developed and some applications 
presented in Refs. 3 and 4. 

By axis translation and successive transforming, other processes 
such as cross-correlation, auto-correlation, convolution, spectral 
analysis, and pattern analysis are possible. Many of these are discussed 
in Refs. 3, 5, and 6. Specifically such systems have found application 
in such diverse areas as character recognition [3 and 7], image improve- 


ment [8], synthetic aperture radar [9], and antenna simulation [6 ]. 


Cwm SGORKeORTHE GHESIS 

This thesis was conceived as a study of the merits of a low cost 
coherent optical processing system for student research. Although most 
systems in use perform two successive transforms, a Single-transform 
system was considered more basic to the problem and such a system was 


constructed. Figure 2 is a photograph of the system. When completed, 
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the system was used as a tool to investigate two otherwise unrelated . 
topics; These topics were a study of Mim ordinvemaracteristics anema 


proposed method of Ship detection via optical detection of the ship wave, 


Dy THESIS OUTLINE 

Section II considers the effects of film in a Single transform optical 
system. A model ofjthe tilm Graimees) proposed and a statistical distm— 
bution of grain size is assumed. From this, a semi-quantitative mathe- 
matical prediction of the rade tansiorm is developed, The grain size 
distribution of Polaroid Type 55P/N film was determined with photomicro- 
graphs and a prediction of the grain transform was made. The grain 
transform prediction was verified with an actual transform of a film 
specimen obtained with the constructed optical system. 

A method of detecting ship waves is considered in section III. A 
relationship between light reflected from a ship wave to the camera and 
the optical transform of the resulting photograph was developed sufficiently 
to verify that a deterministic relationship does exist. A sample ship wave 
photograph was then operated on by a single Fourier transform with matched 
filtering in the frequency plane. The matched filter was rotated to inves- 
tigate the-possibilities of further electrical filtering of the frequency plane 
OUP UL. . 


some suggestions for further research are given in section IV. 
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II. PROPERTIES OF FILM IN A SINGLE TRANSFORM SYSTEM 


A. INTRODUCTION 

Photographic film is the most commonly used medium for the input 
and output in an optical system. The primary advantage of film is its 
high information capacity which Jones [10 ]has compared with other media 
considered as information channels. The three film characteristics 
considered in this section may be thought of as disadvantages or limiting 
factors. The first factor is the time required to process the film, which 
in some cases may run several minutes both for the input and output. 
There are presently two means available to reduce this time; a Polaroid 
type film has a total processing time measured in seconds, while most 
conventional films can be processed in a specially prepared solution 
called a monobath which may require a few minutes. Based solely on 
time, the Polaroid type film seems the obvious choice. Although Polaroid 
negatives require clearing and emulsion hardening before use, the process 
could easily be mechanized and the time made very short. The Polaroid 
Type 55P/N film which provides both a positive and a negative was select- 
ed for evaluation in this project. 

The two other limiting factors are the film's grain structure and its 
non-linear recording characteristics. These two factors were investi- 
gated in particular for Type 55(negative) film, but the method is quite 


general. 


Zl 


B, A FILM GRAIN MODEL 

Photographic films are generally constructed as shown in Figure 3. 
The process of exposing, developing, and fixing the film has the end 
result that the exposed grains are opaque and the unexposed grains have 
been washed from the film. This is described more fully by James and 
Higgins [11]. Even with no exposure, there are a certain minimum 
number of grains which will remain on the film aS opaque spots. This, 
together with the density of the base, sets a minimum achievable density 
which is often referred to as fog level. The number of grains present 
increases as the exposure level increases and achieves a maximum den- 
Sity at some level dependent upon the particular film. Since this grain 
will be Fourier transformed and imaged along with the signal, its effect 
is of some importance. 

The effects of grain have been treated in several papers such as 
Refs. 12, 13, and 14, where the grain is considered as noise having 
statistical properties. The problem is then one of a signal in the presence 
of noise and the techniques of communication theory are applicable. This 
treatment is ideally suited to double transform systems. However, this 
section is concerned with a single transform system, and the spatial 
frequency content of grain is of more interest. 

A review of literature disclosed no appropriate model of film grain 
for a single transform. It became necessary to develop a model and test 
its validity in a single transform system. This model is developed in the 
following paragraphs and though it did not result in a completely mathe- 


matical expression, the result is seen to be of value. 
ZZ 


light sensitive grain 


emulsion 


trans- 
arent 
base 





Figure 3. Cross-section of a negative film. 
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Figure 4 is a photomicrograph of Type 55 (negative) film with no expo- 
sure. The density is thus Sek which is specitied 0.18. The photograga 
suggested a grain structure model of overlapping circular Stops. Since 
the diffraction pattern due to apertures distributed :identically to the 
stops will be essentially the same, the model is initially composed of 
apertures which is mathematically easier to handle, 

Francon [2 | considered the Fourier transform of several apertures 
identical in shape. The model based on Figure 4 is composed of circular 
apertures of varying size but with the size-frequency distribution constant. 
A grain class was defined as a set of grains (or apertures) having one 
specified radius. The work of Francon may now be used to find the trans- 
form of each class. Thus: 

N, - (XaXin + Y2Yin) 

f. (Xz, ya)=f(xz2,ye), L e (2 

n=] 
where N, is the total number of grains in class a, f(Xs,ys), is the Fourier 
transform of one aperture of class a, f,(Xs,ye) is the Fourier transform of 
the set called class a, k equals 27/X, and F is the transforming lens 
focal length. In words, the output due to the set of apertures in class a 
is the summation of N, outputs, each identical in magnitude but of varying 
phase. The total output of a given area of film becomes the summation over 


the area of all grain classes. 


f (2, Ye) = 2 f, (Xo, Yo) (2=2) 
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Figure 4, Photograph of Polaroid Type 55 
(negative) film with no exposure. (1000X). 
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where a is the total number of classes, and (Xe, Ye) is the total Fourier 


transform. The intensity in the transform plane is: 


1! 5, a | £.(x2, ve) |* 
a “a3 
all NPS J FE (2% ntvevin) _ ‘ 
=|E [fe,yeh Ee fe (2-3) 
a=] gts! 
a _k : 
T Ne N, J =| AOS hee Syra) eee (Yin-Yam | 
= 2 \f (Xo,Ye)e L Saad } 
sii n=l m=l 


For N, large. there are a large number of sine and cosine terms, the sum 


of which goes to zero except when m=n. Thus, this term sums to N, and, 
ao | 
I = x | (Xe, Vo). N, | F (2-4) 
© a=l 


But this is not the complete story. From equation (2-3), the intensity at 


the X5,Ye ~Crigiyis: 


ar elle Arp 
1 =r |i? (e-ya)a E > 1 |= . lf (xe-¥aaN,? | We 
p _ oe 2 

a= phe) a= 


The intensity at the origin is a function of the square of the number of 
grains rather than the number of grains as in equation (2-4). Francon [2] 
indicates that, for equal apertures, experiment shows this bright central 
region to be surrounded by a dark .ring. Since there are a large number 
of apertures in a limited surface, there is necessarily a certain distri- 
bution of separation between each grain and its nearest neighbors which 
is not entirely random. The resulting order in the separations produces 
this dark ring. In the same reference, Francon shows that the remainder 


of the output is granular in appearance since the intensity deviation due 


26 


to the phase terms is of the same order as the mean. This effect is noted 


in many other papers and is the cause of the granular appearance which 
is characteristic of coherent illumination. 
The transform of a single circular aperture was then considered. This 


is covered in some detail by Goodman [3] and is: 


f(X2,Yal,e = o* (circa, )] = (ta,)° es f. (2-6) 


where a, is the radius of apertures in class a, J, (2) is the first-order 
Bessel function, and f, is the incident light function (a constant). The 


other terms are defined as follows: 


| =. Sa for Xi oon 
circ(a,) = ie elsewhere 
21 , 
Z=>5 Op (p" = xP+y2) 


(This function has circular symmetry about the optical axis and is sketched 
in cross-section in Figure 5). The intensity distribution due to one 
aperture is the square of the function in Figure 5. Figure 6 is the trans- 
form of a 100 micron radius aperture and was obtained with the experiment- 
al system described later. From equation (2-5), the expression for 


intensity now becomes: 


a 
fh v=) 
I =I x N,ta, e3 (2-7) 
p fe) | Z, 


This neglects the central region and the granularity effects. 
The original intent was to follow this line of analysis to find an 


expression for the total transform with the size-frequency distribution of 


ay, 





Figure 5, The first-order Bessel function. 
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Figure 6, Fourier transform of a 100 micron 
radius aperture. (focal length=50cm), 
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a, Known. However, this involves the distribution function of the out- 
put which is obtained by finding the inverse of the term 2J, (Z)/Z. No 
practical method of obtaining the inverse was found. The output can be 
obtained rather easily, though, for a finite number of classes by graphical 
Summation. This was the technique used to complete the analysis. 

Up to this point, circular apertures have been considered, whereas 
the grain structure model involves circular stops. Defining fs, as the 
transform due to a circular stop model, fs, as the transform due toa 
circular aperture model, and f, as the i che of the incident light 
(source light). the following relationship applies: | 


>» & 


fog (Xp, YE ya)) = tes, ve) (2-8) 


% 


. 


The function f,(Xz,ye) will be the transform of a circular aperture defined 
by the system's limiting aperture. . This is generally the beam size itself, 
Gecince liMicmange Comlpercd LO.a Olaln Size ame pattem twill be 
confined to the region near the origin. At a small distance from the 
origin, only fs,.and fo, are Significant and: 

foc (Xo,Yo) = - fog (K2g,Ye) 8) 
Ol: 

ined (Xo,Y2) = fos (X2, Ya) (2= iG) 
That is to say that the intensity distribution of the transforms of comple- 
mentary inputs are identical except at the origin. This is Babinet's 
Complementary Aperture Theorem [2]. 

Considering all of the above factors, the film grain transform will 


exhibit a brilliant spot at the origin due to the function f, (xs,y.) (in 
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effect - the DC term), and the N,* factor of equation (2-5). Surrounding 
this is a dark annular ring as noted above. The remainder of the pattern 
will be granular in appearance with a mean distribution given by equation 
iZ-7). 

Thus far this analysis has ignored the fact that the film thickness 
varies. Although some variations will occur solely because of manu- 
facturing tolerances, the most prominent effect is due to the processing 
action of exposed grain removal [3]. By defining 0 (X14 + Yan) to be the 


additional phase introduced at each grain center, equation (2-1) becomes: 


Tk 
Na a (X3X1 nt+Ve2VYin) a bn | 
£.(x2,Va)=! (Xo, Vole Ji e (2-11) 

n=1 
The analysis then follows as before with a similar result. The output 
pattern is essentially the same as given in equations (2-4), (2-5), and 
(2-7) except near the origin. Notethat for x,=y.,=0, d is not necessarily 
equal to zero, and some of the intensity at the origin must be dispersed 
over a Slightly larger area. The dark ring surrounding the origin will not 


necessarily be at the same radius since the causal order depends on both 


the grain-center separations and the emulsion thickness variations. 


C. THE PREDICTED TRANSFORM RELATED TO FILM DENSITY 

A size-frequency distribution was obtained from Figure 4 and the 
results given in Table I and Figure 7. A size count such as this is some- 
what arbitrary and very dependent on the judgement of the observer. The 


guidelines used were: (1) for clumps of more than one grain, choose the 
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Figure 7, Size-frequency curve of 
Polaroid Type 55 (negative) grain. 
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TABLE I, GRAIN SIZE COUNT. 


Class Diameter 
Diameter Range Number Weight (number/total) 
0.4} 0.3-0.5p Zeg 234 
Oi 0.6-0.8 Bie .341 
luge 0,9-1.1 Salis 2:85 
ie 1. Zee 84 .076 
Nels. lL. .oF 1.7 47 .042 
1.9 eee2 .0 les: .014 
kale Ca 3 4 .004 
Cas 2.4-2.6 = Z000 
Cae Ot 2./- @ 0 .000 
Total 1106 1 #001 


largest Sub-units which could be considered as circular and count each 
as one grain, (2) whenever possible, choose sub-units such that overlap 
is minimized. The resulting curve of Figure 7 is Similar to a typical 
curve given in Ref. ll. 

The data of Table I was used to determine the intensity distribution 
in accordance with equation (2-7). Figure 8 is a plot of the transforms of 
each size-class and the total transform. There are no sharply defined 
nulls or peaks in the total transform which was to be expected since this 
is the transform of a continuously distributed input. The curve is a 
smoothly decreasing function of radius, which is difficult to study by 
experiment without some means of accurate light intensity measurement. 
Photographic film does not supply the desired accuracy and, at this time, 
no other means was available. At this point, it seems worth mentioning 
that the transforms of many signals of interest are not continuous Spectra. 


As such their patterns will approach point spectra which are more easily 
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Figure 8, Fourier transform of Type 55 (negative) film grain 
by graphical summation. (focal length=10 cm). 
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analyzed. Film can provide only a rough quantitative analysis for the 
continuous distributions as in the present case. 

Type 55 (negative) film was chosen for the output medium. The 
D/logE curve (see Appendix A) is non-linear except in the central region 
of the slope. This is generally true of any film. In the linear region, 


the following relation holds: 


D = y log E + constant 


y log IT + constant (2-12) 
where E is the photographic exposure, and T is the time duration of the 
exposure, This relation is valid from D=0.375 to 1.5. The range of 
output intensities far exceeded the linear range of the film and many 
exposures were needed at various values of T to cover the range. Thus, 
by considering the linear density regions of output films exposed at 
various values of T, the effective linear density range was extended. 
Since the intensity at a given radius must be a constant for all photo- 


graphs, two successive photographs were compared by letting: 





= I, Ty 
AD = y log LT. 
aMeesct I,=[,. Thus: 
qT, 
AD = » log.—— (2-13) 
Tp 


In practice, a AD was determined by comparing densities at equal radii 


in two photographs. 
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D. THE GRAIN TRANSFORM BY EXPERIMENT 

The optical system for this experiment is diagrammed in Figure 9, 
The position of L, and P, are reversed compared to the normal arrangement, 
which alters the transform kernal from exp(jk/F) (x,xX2gty, ye) to exp(jk/d,) 
(x,X=ty,ye) and offers some control over the size of the output [3]. The 
input no-exposure tranSparency was placed against an optical flat with 
a matching fluid (mineral oil) between the two. Several exposures of the 
transform were taken to cover the intensity range with the input trans- 
parency base side against the fluid. A similar set of exposures was made 
with the emulsion side against the fluid. A Single exposure was taken 
with the input transparency unmatched by the fluid. A visual comparison 
was made with the single unmatched output and the two matched cases at 
the same exposure level. The case with the fluid against the film base 
appeared virtually identical in density to the unmatched case, indicating 
that the base thickness variations have a negligible effect for this film. 
The emulsion-matched set of exposures correspond to equation (2-7) with 
no random phase variations. The base~matched set of exposures corre- 
spond to equation (2-11) because of the phase variations due to the 
unmatched emulsion side. 

The density of each of the exposures was recorded uSing a micro- 
densitometer to track a diametric line through the origin. Figures l0and 
11 are typical density recordings from each of the sets. Figures 12 and 


13 are the corresponding Type 55 positives for these two recordings. 
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Figure 9, Optical system for grain structure Fourier transform. 
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Figure 13. Transform of 
grain in Type 55 (negative) 
film with base surface 
fluid matched. (Density 
recording is Figure 11), 
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Figure 12, Transform of 
grain in Type 55 (negative) 
film with emulsion surface 
fluid matched. (Density 
recording is Figure 10). 





The exposure level is not the same for these two photographs and no 
comparison of the exact density levels should be made. They are pre- 
sented here only as typical of each set. In both figures, the peak is in 
saturation while the minimum level is below Dan’ The bright central 
peak and Wark annular ring are obvious in Figures 11 and 13. In Figure 
ica faint ring is visible at a larger radius and appears as two Spikes 
in Figure 1 0. Although this ring may also be due to the non-random 
order of apertures separated by short distances, it would be impossible 
to verify it without further knowledge of the film construction. A fine 
structure is evident in the recordings but the true character of granularity 
cannot by recorded by the microdensitometer since the variations occur 
in the order of one wavelength (0.6328 microns). The "tailing-off" in 
both figures is due to the film non-linearity in the toe of the D/logE 
curve. 

The series of exposures in each set were compared indensity as 
noted above. The resulting density curves are plotted in Figure 14. The 
Eoretical curve of Figure 8 was converted to density levels and also 
plotted in Figure 14. Note that the theoretical curve was derived from 
equation (2-7) which does not consider the effects near the origin. It 
must alSo be in error due to the effect of overlapping grains. When the 
grains overlap, the resulting shape is no longer earentar and the dimension 
in one direction becomes much larger on the order of multiples of the 
smaller dimensions. Referring to Figure 8, note that the transforms of 


the larger grains affect the resulting curves at low values of radius. The 
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Figure 14, Normalized Fourier transform curves for Type 55 

(negative) film grain. Density levels based on a linear 


D/logE characteristic using Type 55 (negative) to record 
the output. (focal length=10cm). 


42 


small grain transforms have a proportionately greater effect at large 
radii. The photomicrograph of Figure 4 shows very little overlapping for 
the smaller grains, keeping in mind the judgement criteria upon which 
the grain count was based. Thus, the experimental results at large 
radii should closely approximate theory. The curves of Figure 14 were 
normalized by equating the density levels ata radius of 12 cm. 

The experimental and theoretical curves show excellent agreement 
between 5 and 12 cm radius. There is no reason to believe that they 
would not continue to agree at larger radii but the curve could not be 
extended further because of lens curvature of field. The density departed 
from the theoretical curve for small radii in the manner predicted. That 
is, the intensity was higher near the origin due to the strong DC term 
and the effects of overlapping. It was noted also that the peak of the 
base-matched curve was less than for the emulsion-matched curve and 
that the intensity was somewhat more dispersed about the origin. 

The method of this section has been both quantitative and qualitative. 
Exact agreement between theory and experiment is not possible because 
the model is not an exact copy of the true grain structure. Had a better 
means of determining the intensity distribution been available, the 
normalizing procedures would not have been necessary and the results 
would, admittedly, carry more validity. However, the close agreement 
which does exist between theory and experiment is sufficiently convincing 


to prove the essential correctness of the model. 
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E, THE SIGNAL-PLUS-GRAIN TRANSPARENCY 

Thus far only the transform of the minimum grain level has been 
considered. This must now be related to a transparency composed of 
this minimum grain density with a signal superimposed on it... A,typiteal 
transparency may consist of a continuous distribution of densities ora 
set of binary density levels, i and Dan! but in either case the 
Signal is on a background field of Ces or Da = Since there exists a 
limited density range, the signal itself must consist of Some grain 
structure even in the binary case. 

Figures 15 and 16 are photomicrographs of Type 55 (negative) film 
with density levels greater than Di and they should be compared with 
Figure 4, It is obvious from Figure 14 that the circular grain model is far 
less representative of the actual grain structure when the grain structure 
is composed of irregular clumps. The nature of the exposure-development 
process is not fully understood, though several plausible theories have 
been proposed [11]. It is important here to note only that a clumping 
effect does occur and that it increases with density. Larger dimensions 
are now involved and, from the previous discussion in the effect of 
increased dimensions, an increase in intensity near the origin occurs at 
the cost of intensity at large radii. Assuming the input transparency 
area to be large compared to the size of any one grain, the size distri- 
bution is identical in all directions, and the output transform due to grain 
will still be symmetrical about the origin. The central DC term is reduced 


at the same time due to the overall increase in density. 
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Figure 15, Photomicrograph of Type 55 
(negative) film. Density 0.7. 


Figure 16. Photomicrograph of Type 95 
(negative) film. Density 1.7. 
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The shifting of energy to smaller radii due to.larger grain dimensions 
continues until the total grain area of the negative is dp plot taal equal 
to the transparent area. Further increases in exposure produce negatives 
whose grain structure iS approximately the Fomolenen t of that for the 
negatives with lower densities. Thatis, the structure is one of trans- 
parent apertures rather than opacities stops. Indeed, in Figure 16, the 
structure appears to be a rough Pomplement of Figure 4. The Complement- 
ary Aperture Theorem of Babinet predicts identical diffraction patterns for 
complementary Screens except for the central DC term, Exact complements 
oe be expected because the apertures of the high density negatives 
do not have the circular structure of the stops of the low density nega- 
tives, ey NCEE structures, an increase in density results in a shift of 
energy in the transform from small to large radii. The central DC term 
continues to decrease, 

A signal transparency will contain more than one density level and 
the output transform due to the grain alone will be dependent on the 
existing input density levels and relative areas which they occupy. The 
binary input consists of areas of complementary grain structures and the 
output will be the signal transform ate a grain transform which approaches 
that in Figure 8 or 12. The worst case occurs when the Wout consists of 
intermediate densities where the grain size is about equal to the aperture 
(between grain) size, The transform intensities are more dispersed. 

Returning to Figure 14 and remembering that photographic film has 


a limited recording range (AD is about 1.5 for Type 55), it is seen that an 


adjustment of exposure level to place the tail of the curve below Din 


4 


leaves only the grain transform near the origin. The effect of grain (with 
the emulsion matched) in the transform would be a central disc on the 

order of 0.2 cm radius for a focal length of 10 cm or, in general, 0.02X 
(focal length in cm) for Type 55 (negative) film. This means that the signal 
can have no significant frequency content less than about 2x10° cycles/cm 
for this film. 

Figures 17, 18, and 19 are a series of signal transforms which were 
designed to test the above hypothesis and serve as a visual demonstration 
of transforms in the presence of grain. In practice, the best output was 
obtained at an exposure level which included some of the tail of the 
density curve of Figure 14. Figure 17 is the transform of Din grain 
alone. The signal input for Figures 18 and 19 were 2mm high characters 
with density Din against a field of Dax’ The fine structure of con- 
centric rings seen in all three photographs is due to a small aperture 
mask required to select a particular character on the input transparency. 
Although the grain transform does degrade the image, the distinguishing 


features of the Signal transforms are readily apparent. 


PF, CONCLUSIONS 
Although the experimental grain transform differed significantly 
from the theoretical grain transform, the essential correctness of the 
grain model was verified. The overall effect of the grain appeared as a 
reduction in bandwidth by setting a lower limit of useful frequency response 


in cycles/cm. This effect appears to contradict intuition since grain is 
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Figure 17. Trancnor meen 
film grain. Exposure is 
the same as for Figures 
18 and 19. (focal length 
is 50 cm). 





Figure 18, Transform of 
the character '4' in the 
presence of film grain. 
(focal length is 50 cm). 
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Figure 19, Transform of the character '5' 
in the presence of film grain. (focal length 
eS Snel oleanh 
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usually thought of in terms of limiting the fine detail, or the high Spatial 
frequencies in an image. However, the following discussion shows that 
this 1S i special case and the results are somewhat misleading beceuse 
a relatively grainy filmwas used. | 

It must be noted that no part of the signal, including the fine detail, 
can have dimensions less than the diameter of the smailest grain. In 
most cases, the ewest spatial frequency of even the finest Retail in the 
signal will fall in the central region where the grain transform is Signi- 
SOG Jai bi The higher frequencies of the signal are visible here only because 
the signal transform intensity is significantly greater than the grain trans— 
form intensity in this frequency range. As seen earlier, the total inten- 
sity Sime grain Geach is a function of the total number of grains (or 
apertures in the case of 1D) ae in the input transparency. The grain 
transform intensity was reduced by limiting the area of the input trans- | 
parency to include little more than the signal area as in Figures 17, 18, 
and 19. Retention of the low frequency signal spectrum, including that 
in the fine detail would require that the signal transform intensity be 
greater than the transform intensity over the range of the Signal bandwidth. 
Extremely fine grain films are required to disperse the grain transform 
spectrum and "flatten" the curve of Figure Ge 

Consider a limiting case, shown in Figure 20, in which the grain 
transform is so dispersed that it appears essentially constant over the 
spatial frequency bandwidth of the signal. Ina typical signal, the 


transform intensity would be a decreasing function of frequency. Above 
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Figure 20. Idealized signal and grain spectrums 
for a signal on fine grain film, 
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Figure 21. Idealized signal and grain spectrums 
for a signal on coarse grain film. 
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some high frequency limit, the Signal spectrum intensity will be less than 
the grain spectrum intensity. Under these conditions, the grain structure 
is seen to have a greater effect on the high frequencies or the fine detail 
of the Signal. Figure 21 illustrates a case sales to that of Figure 20, 
but one in which the grains are coarse. This case resembles an experi- 
mental — . In Figure 21, the grain spectrum intensity is greater 

than the signal spectrum intensity at both low and high frequencies and 
effectively creates a pass=band, The high frequency limit was not ob-wi\ 
served in Figures 18 and 19 because the intensity range in the output 
exceeded the film density recording range. Further, since the grain 
transform intensity obscured most of the important lower signal frequen- 
cies in Figures 18 and 19, reconstruction of the signal by Fourier trans- 
form of the output transparency was found to be impossible. 

Thus, film grain in a coherent optical system has an effect similar to 
noise in electrical circuits, The film grain spectrum in the limiting case 
of Figure 20 appears similar to GausSian noise in that its spectrum is flat 
over the range of interest, It appears that in many cases involving 
binary signals on fine grain film, this noise could be assumed Gaussian, 
at least aS an approximation. In the discussion relating to Figures 15 
and 16, it was noted that the film grain spectrum is a function of the film 
density for densities other than aren and Dax’ This would mean that 
there would necessarily be some correlation between grain noise and a 
Signal comprising intermediate densities. The assumption of Gaussian 


noise could then be, at best, an approximation. 
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ae INTRODUCTION 

The wave pattern produced by a pressure disturbance moving on 
the surface of the water was first considered by Lord Kelvin [15] in 1891. 
For a disturbance such as a Ship moving ona Straight course, the wave 
height is considerable inside an angle bounded by two horizontal rays at 
419.5 degrees with respect to the ship's course line [16]. Further, the 
angle is not a function of Speed. Such a pattern, when considered 
against a random sea wave pattern, is unique and optical Fourier trans- 
forming techniques may be applied for its detection. 

The method used in this project was the single transform of a wake 
photograph with a matched filter in the first transform plane. The matched 
filter was rotated on the optical axis and the resulting transmitted light 
measured with a CdS cell to examine the possibility of further electrical 
filtering. A thorough investigation was not conducted since this was 
considered as a feasibility study and as a practical use-test of the optical 


system. 


B. THE WAVE PATTERN 

The free surface equation of ship waves is very difficult to derive 
when exact hydrodynamical theory is used. The solution of Lord Kelvin's 
as well as later work are approximations [16 and 17]. Even then the 


results are quite involved and not completely germane to the problem of 


Bye) 


detecting the waves. Ursell [16] has predicted and plotted contours for 
a ship’s wave near the critical angle, @=19.5 degrees. Figure 22 is one 
of these plots. 

Along the critical angle, the wave crests #. inclined at an angle 
of 54.75 degrees to the track of the ship. Outside the critical angle, 
the wave height falls rapidly to zero. The horiZontal angle of inclination 
between wave crests and the track decreases from 54.75 degrees to zero 
degrees for angles from less than 6, to points on the track, However, 
the wave height also decreases inside 8, although not as rapidly as out- 
side 9@,. Thus there is some horizontal inclination angle between 54.75 
degrees and zero degrees at which the wave heights become insignificant. 
By way of comparison, the wave heights near the critical angle die out as 
an inverse function of the cube root of distance while, in the interior, 
they die out as an inverse function of the square root of distance [18]. 
The wave height and period are functions of the ship's speed and, from 
Figure 22, it can be seen that the period also decreases with angle. In 
Summary, as the angle decreases from 9, to zero, the wave height 
decreases, the wave period decreases, and the wave crests' angle of 
inclination with respect to the track decreases (normal angle increases). 
These factors allow a qualitative prediction of the Fourier transform of 
such a pattern to be of the form in Figure 23. oN is the direction of wave 
travel (normal angle) for wave components near @,, and * is the direction 
of wave travel for wave components which exist near the track. oa equals 


90-54.75 degrees or 35.25 degrees. 
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Figure 22, Contours of equal magnitude of the dimensionless 
surface elevation pV C., (x,y)/2gP. due to a concentrated force 
Py moving with velocity V, for gr /Vo=2 0 r/d near 100. (Here 

r denotes the distance from the disturbance, the wavelength 
of waves moving with phase velocity V.) For the sake of 
clarity only positive values of the elevation are shown, [After 
Ursell, Ref. 16]. 
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Figure 23, Approximate Fourier transform 
of a ship's wave, 
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C. PHOTOGRAPHIC PROPERTIES 

Stillwell [19 ] analyzed the relation between the light scattered from 
asea surface and the optical density of a scene photograph. The method 
is outlined here with some changes which are ——e to the present 
problem. 

Let the ship's wave height pattern be described by H\ x,y), where x 
and y are taken to be on the sea surface, Assume the sea surface is 
illuminated by a completely clear or a completely overcast sky. The 
illumination at the camera becomes: 

M(n) = KN(y) P'(B,) (3-1) 
where K is a constant, L(@ the luminance of the sky in the direction 
reflected into the camera, and [(§8,) is the reflectivity at the incident 
angle 6, . Figure 24 describes the geometry for a one-dimensional surface. 
For the two-dimensional surface, 8 and @ are simple functions of 6 and 
O@ . The water reflectivity is a function of the incident angle and this is 
illustrated in Figure 25. Thus the camera illumination under these condi- 
tions is a function of the sea surface wave height structure, H(x,y). The 
camera illumination can now be described as: 

M(x, y) = K,H{x,y) (3-2) 
where K, is a constant. (Note: The illumination is also a function of 
the angle over the field covered by the photograph. The illumination is 
maximum directly away from the observer and is present in all sea photo- 
graphs [19], A narrow field of view will be assumed so that this factor 
is eliminated.) Using equation (2-12), the recorded film density is: 


D(x,,y,) = ylog I(x,,y,)T + Kz (constant) (3-3) 
96 
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Figure 24, One-dimensional sea 
Surface reflection geometry. 
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Figure 25. Water reflectivity. In the visible 
region, water reflectivity is related to the 
angle from the normal by the above curve, In 
addition to the wavelength dependence of the 
reflectivity, the vertically (V) and horizontally 
(H) polarized components have different reflec- 
tivities. [After Stillwell, Ref. 19]. 
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where T is the exposure time and the coordinates , Cys are in the 
film plane. Figure 26 illustrates the relation between the surface co- 
ordinates, (x,y), and the film coordinates, (x, ,y, ). The magnification, 
figs a funetron- or tic illimetorob)ecmal Stance _— r,t+r, is 
approximately equal to r,+rz over the region recorded, then m may be 
assumed constant over the region and simple trigonometric relations may 
be applied. This approximation is generally true since the camera height 
is usually much greater than the focal length, The coordinates, x and X,, 
are arbitrarily assumed to be in horizontal planes such that the relation, 
x, —Mnxyeas , SiMirercoordimeloos y andayware thenin a vertical plane 
containing: bethesy andsy,. The plane, P', is drawn nommalhto therlinese! 


Sight and the coordinates related to this line become: 


and 0 <m <l (3-4) 
The following angles may be derived. 
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And finally, the Law of Sines, 


y,/m —_p pee (3-5) 
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Figure 26, Sea surface-film plane geometry. 
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gives the required relation. 


iy oo (3-6) 
1 + We 


For a high altitude camera Such as ina satellite, g is approximately 90 


Vuge= 


degrees and equation (3-5) becomes, 


Vemma) Y 


and (3-7) 


D. THE SHIP WAVE TRANSFORM 

The exposed photograph contains the sea surface information in the 
form of density variations as given by equation (3-3). This transparency 
is inserted into the input plane of a single Fourier transform system. The 
transmitted light amplitude is given by: 


ai. a 
Dey ’ Yi) 


fixwy, ) eR, ee (3-8) 


The light intensity function in :the output or transform plane is; 
I erve) = ke FEF, Oy.) 1! (3-9) 


Combining equations (3-2, 3, 8, and 9), the final expression is: 


1 %arVo) ae i LHx,y)7 0/24 - (3-10) 


where: K, = Ka(f, 10°?K,T)? and x,y are 
given by equation (3-6) or (3-7) 
Although the above method is not complete and includes several assump- 
tions, the results are sufficient to verify the existence of a predictable 
relationship between a ship wave and the Fourier transform of its .. ~ 


photograph. 60 


In a practical Situation the ship wave must be considered in the 
presence of sea waves. The sea waves will not have the unique Fourier 
transform of the ship wave and, on the basis of this alone, should be 
easily separated in the transform plane. However, the sea waves do 
contribute to an overall noise level and it may be assumed that there is 
some minimum ratio of Ship wave height to sea wave height below which 
the optical detection would be useless. Later mention will be made of 


this in section IV, 


E. ASHIP WAVE DETECTION EXPERIMENT 

The photograph of Figure 28 is a ship wave ina low sea state. The 
characteristic features are readily apparent along the critical angle. The 
transparency for Figure 28 was inserted in the input plane of the optical 
system illustrated in Figure 27, and an exposure was made in the transform 
plane. This is Shown in Figure 29. The transform was similar to that 
predicted in section IJI-B. The differing density of the two pairs of wave 
lines may be traced to Figure 28 where the location of the sun produced 
differing modulation levels on the two sides of the wake. It was noted 
that this affected primarily the density rather than the Shape of the out- 
put transform. The remainder of the pattern is composed of the grain 
transform and the central DC intensity considered in section JI, and the 
transform of the sea waves which appear as noted in Ref. 19. 

A crude matched filter was produced by cutting out the portion of 
Figure 29 which represents the ship wave alone. This filter was then 


aligned in a mount which allowed rotation of the filter about a plane 
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Figure 28. Ship wave in 
a low sea state, 





Figure 29. Fourier transform 
of the ship wave of Fiqure 28, 
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normal to the optical axis. The mounted filter was placed at plane, P2, 
and a Cd& cell covered with diffusing glass placed behind it. The geom- 
etry was designed so that the area of the transform pattern matched the CdS 
cell face. The cell output was then readasa EeRetion of angular filter 
position. 

This particular experimental arrangement was chosen to investigate 
the possibility of performing additional filtering by rotation of the matched 
filter at a constant rate. The CdS cell output would then be a time cor- 
relation of the signal and the matched filter. However, it must be noted 
that a CdS cell has a slow response time and could not be used at prac- 
tical rotation speeds. Referring to Figure 29, the output should consist 
of an average level plus components due to Sea waves and the ship wave 
with a fundamental frequency of twice the rotation rate. The ship wave 
component would have a shorter time duration and produce different fre- 
quency harmonic terms. The resulting data is presented in Figure 30. 
This is the best run of a series. The data not included here was similar 
to that of Figure 30 but, in general, was degraded by inaccurate filter 
construction and filter centering. The construction problem could have 
been eased by uSing a photographic transparency filter, however, this 
would also have introduced undesirable film grain effects. The portion 
of the curve in Figure 30 representing the ship wave is quite prominent 
and suggests the possibility of visual detection. However, since the 
intended use of the output would involve frequency analysis, the curve 


was analyzed further for frequency content. 
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The FORTRAN program, FORIT, which computes Fourier series coef- 
ficients for a tabulated, periodic function was used for frequency analysis. 
Nineteen equally spaced points over 180 degrees of Figure 30 were used 
to obtain coefficients to the ninth harmonic, For comparison a similar 
curve, in which the ship wave component was replaced by a straight 
line joining the component end-points. The resulting exponential Fourier 
series coefficients are presented in Table II. The values normalized to 
the DCoclevel 


TABLE II, EXPONENTIAL FOURIER SERIES 
COMMA MAINS Clee abla aa @xelwe Glo LIL (Owar ae. 


HARMONIC COEFF. WITH COEFF. W/O NORMALIZED NORMALIZED 


WAKE WAKE WITH WAKE W/O WAKE 
0 Z onal DOS SEV 1.0000 Levooe 
] URO9Z7 0.2084 Ossi 0.0927 
2 me ov 00833 0.0373 0.40397 
3 eli O59 0.0881 O20 13 0. @3gZ 
4 Vee 5 Ooo 0.0426 0.0307 
S 0.0726 0.0263 0.03 a O au 7 
6 0.0679 OnmeGsy'Z O@20d 0.0174 
7 0.0642 0.0486 0.0275 0 .@2Z16 
8 0.0490 0.0038 OS0Z2 20 0 G0 lg 
3 0.0543 OU 178 UZ 3 0.0080 


are given in Table II and plotted in Figure 31. A marked difference was 
noted particularly in the fundamental and the second and third harmonics. 
It was concluded that the output is Ssusceptable to frequency filtering to 


indicate the presence of a ship wave in the input transparency. 


FP S=GONGBUSTONS 
Thus, the combined optical-electrical detection of a ship wave has 
been shown to be feasible, at least on a crude level. Since the transform 


of a ship's wake is always centered on the optical axis, the system is 
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not dependent on the wave position in the input field. The rotation of the 
matched filter eliminates the dependence on angular positioning. A phase 
reference derived from the filter rotation drive could conceivably be used 
to determine the ship's course with a 180 degree ambiguity by comparison 
of the phase of the ship wave component in the output signal to the phase 
reference. These factors suggest the use of such a system in an auto- 
matic alarm system to indicate the presence of a moving vessel in some 


monitored area of the ocean, 


IV. SUGGESTIONS FOR FURTHER RESEARCH 


A, FILM CHARACTERISTICS 

The film grain model developed in section II allowed a semi-quanti- 
tative prediction of the grain transform. The ultimate limitation of this 
method lies in the grain model itself. However, the results would have 
been more useful with a closed mathematical expression for the total 
transform with the input size-frequency distribution known. As stated 
before, from equation (2-7), this involves the distribution function of the 
output which is obtained by finding the inverse of the term 2J,(Z)/Z. The 
expression for the grain transform could be obtained in a closed form if 
this stumbling block were overcome. Since the input size-frequency 
distribution is often known for a type of film, the results could be applied 


almost directly to predict the effects of grain for any film type. 


B. SHIP WAVES 

The ship wave detection problem has been treated here on a fairly 
basic level with several approximations and assumptions. Since the 
concept does have obvious applications in surveillance of the seas, it 
seems worthwhile to define the problem more fully. Some of the problems 
involved are: (1) the variation of reflected light modulation with the sun 
and observer positions, (2) the relation between a ship wave and sea 
waves travelling at various relative surface angles, and (3) the signal 


to noise ratio, 
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C. SEA SURFACE ENERGY SPECTRA 

In Ref, 19, Stillwell describes the use of a Single optical Fourier 
transform system to determine the energy Spectrum of sea surface wave 
heights or wave slopes. Research into this area is cof importance not 
only to oceanographers but is also vital to studies in radar scattering and 
a@eustic propagation. 

A brief outline of Ref. 19 is given here. As noted in section III, the 
light arriving at a camera above a sea Surface is a function of the wave 
slopes because of the incidence angle dependence of water reflectivity. 
The density of the exposed film is then determined by the sea surface 
characteristics and the problem geometry. The optical Fourier transform 
of such an exposed film can be expressed in terms of the film density 
and source light intensity. Since the transform output records the spatial 
frequency components of the input, it will also represent the sea surface 
spatial frequency components with a scale or magnitude factor involved. 
It is shown that the spatial frequency components of the sea surface are 


related to the slope energy Spectrum, which is expressed as: 


&(K,) = C 10 Y sec? yb (4-1) 


where C is a constant, D(k/m, w) is the output or transform photograph 
density function, k is the radian wave number, mis the magnification 
meleror, t) is an angle between the obServer=suriace normal plane anemic 
observer-zenith plane. w= 0 ona line directly away from the observer. 
Stillwell presents Some results by plotting density contours of the trans- 


form output. Atypical contour is Shown in Figure 32. 


69 


WAVELENGTH [cm] 


a 
oO 
oO 


> Oo in 
Oo oO o 


S2 





Figure 32. Density contours showing the long wave- 
length spectrum and spectral peak for a light to moderate 
wind of long fetch. [After Stillwell, Ref. 19]. 
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peveral attempts at analyzing sea surface photographs were made 
using the present optical system. The original photographs were 8X10 
inch transparencies, of which an area of 45X4¢ inches was selected and 
reduced to 1/20th size, The area selected is shown in Figure 33 and the 
transform output in Figure 34. A microdensitometer was used to obtain 
the contours in Figure 35>) The results were genendliyeunouee es ful but 
not totally so. The contours do indicate some pattern which could only 
be attributed to the sea surface structure, and they are roughly similar 
to those in Figure 32. Most of the difficulties can be traced to the 
quality of the input. The original 8X10 inch transparency was grainy and 
the reduction of a photograph usually entails some loss of high frequency 
content. Both of these factors are apparent in Figures 34 and 35, and they 
are serious enough to make the reduced negative virtually useless. 

It is suggested that further work might be done in this area using 
sea surface photographs taken for the prime purpose of use as inputs to 


the optical Fourier transform system, 
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APPENDIX A, FILM CHARACTERISTICS 
The typical D/logE curve has the form of Figure 36. A standard 
approximation for the negative type film is shown in Figure 37, The result- 
ing equation which describes the curve of Figure 37 is: 


— i | GPA < 
D y [log E - log i] for E(D A) _—< E(D ax! 


sage : 
ees ee : (Dax! (A-1) 


D.., fon L = § eG) 
min min 


Often the spectral sensitivity curve, gamma, D , and D_._ are the 
max min 


only values given. Equation (A-1) must then be used to determine the 
density as a function of exposure. Spectral sensitivity curves relate 
radiant light intensity (or energy density) as a function of wavelength to 
reach a specified film density, usually 0.5. The exposure, E, is the 
incident luminous flux, U, times the duration of exposure. Luminous 
flux and radiant flux, R, are related by [20]: 

U 


luminous efficiency (7) = - (A~2) 


which is 150 lumens/watt for 6328A. The exposure is then: 


® 


p= —t— --2ht_ (A-3) 
unit area unit area 


The constant, i, is evaluated using the radiant light intensity for D=0.5 
from the spectral response curve. The result is; 
D =¥y [logE - log E(D=0.5)] + 0.5 (A-4) 
[Ore Cer <E(D_ ) 


min OX 
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Figure 36. Typical film characteristic curves. 
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1 logE 


Figure 37, Approximation of the film 
characteristic curve, 


When exact light levels are not known or desired, the following form is 
useful. 
D = y log E + constant . (A=5) 

The effect of reciprocity failure is not considered here, Reciprocity 
failure occurs when exposures involve a high intensity and a very short 
exposure time or a low intensity and very long exposure time. The effect 
is a reduction of film sensitivity and the exposure time must be increased 


accordingly [11]. 
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APPENDIX B. LIGHT INTENSITY METER 

The light intensity meter described here was found useful in con- 
junction with the optical system, It can provide a rough measure of the 
intensity distribution at any point in the system and eliminate some ot the 
guess work to obtain a proper exposure, 

The circuit is diagrammed in Figure 38 and is basically a high gain 
operational amplifier with a cadmium sulfide photocell in the feedback 
loop (23 |. This configures was used to reduce the effects of temper- 
= and power supply | tee ‘Transistors Q, and Qs, act as a 
differential amplifier with a my ole ayer Transistor Q, is a DC level 
shifter and amplifier while transistor Q, drives the meter. 

The photocell was mountedMmimga teleseopre wuee behind a 0.5mm 
diameter pinhole. The telescopic tube was affixed to a bracket with 
vertical and horizontal a. so that the photocell position could be 
adjusted over the face of the bracket. The bracket slides into the film 
pack holder with the pinhole in the film plane. The completed unit is 
Shown in Figure 39% 

The unit was calibrated from zero to 2990 meter-candles using 
diffuse illumination and a Weston [llumination Meter. There was no 
standard available to calibrate the meter above 2990 meter-candles. The 
calibration above 2990 meter-candles was accomplished by comparing the 
readings at two distances from the source. Letting I, (d) define a calibrated 


reading at distance d, the relation is: 
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R1-500 ohm pot 
R2-2500 ohm pot 
R3-25,000 ohm pot 
R4-100,000 ohm pot 
R5-10,000 ohm, =W 
R6-3000 ohm, =W 
R7-1000 ohm, +W 
R8, Rl1- 470 ohm, = W 
R9-10,000 ohm, = W 
R10-2700 ohm, = W 
nee CLS0S-L 


Cl-10uf, 15V elec, 

S1-DPST, normally open 
52Z-SP, 5 pos. selector switch 
Bl=-1.5V 

B2-9V 

Q1, Q2-2N3391 

Q3-2N4126 

Q4-2N3391 


Figure 38, Light Intensity Meter circuit. 
[After Barnum, Ref. 21], 
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Figure 39. Light Intensity Meter 
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= constant, K (B-1) 


The distances, d, and dz, were used for all readings and the readings 
outside the calibrated range become: 


ac alienated (d,) = KI. (de) (B-2) 


The resulting accuracy cannot be as close as +5%. The calibration curves 
are givenin Figure 40. A fourth meter scale was not calibrated because 

a light source of such high intensity was not available. This scaleis 
beyond the range of intensities in the optical transform system and was 
not needed. 

In use, the probe is positioned to monitor the light intensity ina 
region of the output bs where a linear density response is desired, An 
approximately correct first exposure can be obtained by exposing fora 
density of 1.0 for most films [11]. Equation (A-4), solved for E to reach 


a density of 1.0, becomes: 


Ee =e {D=0), 5) 19°:5/¥ (B-3) 
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Figure 40. Light Intensity Meter calibration curves. Scale 
readings are: Xl for Scale A(#5%), X10 for Scale B(+5% for 
meter reading less than 7 and 10% above 7), and X]00 for 
Scale C (20%). Functional tests: BATTERY-13.3 meter 
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